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ABSTRACT 

We develop a phenomenological theory that aims to account for the origin of 
the large eccentricities of extrasolar planets and that of the small eccentricities 
in the solar system, the preference for apsidal alignment in non-resonant multi¬ 
planet systems, and the origin of the similarities in the eccentricity distribution 
of extra-solar planets and that of spectroscopic binary stars. We show that if a 
physical process is weakly dependent on the local dynamics of the companion and 
imparts a small relative acceleration to the star-companion system, the eccentric¬ 
ity of the companion’s orbit is excited to large values depending on the direction 
and duration of acceleration. A natural candidate for such processes are asym¬ 
metric stellar jets and star-disk winds. When the acceleration originates from 
a precessing jet, large eccentricities can be excited by the resonance of the jet’s 
precession frequency with the induced acceleration’s excitation frequency even 
for nearly perpendicular jets. Precession also reduces the eccentricity amplitude 
far inside the resonance radius. The acceleration’s strength is best constrained 
in multiplanet systems because of the companions’ mutual gravitational pertur¬ 
bations, while the acceleration’s duration is bounded by the condition that the 
residual velocity imparted to the star is smaller than the stellar velocity disper¬ 
sion in the Galaxy. In the outer parts of the star-companion system where the 
acceleration excitation time is comparable to or smaller than the orbital period, 
significant radial migration takes place which may have important consequences 
for the dynamics of the minor body populations in the solar system. The theory 
is illustrated with the v Andromedae binary system. 

Subject headings: celestial mechanics - planetary systems — stars: mass loss - 
stars: winds, outflows 


1. Introduction 

Since the discovery of the first Jupiter-like planet around 51 Pegasi in 1995, the number 
of extrasolar planets has exceeded 150, some of which reside in 14 multiplanet systems 


- 2 - 


(Schneider 1996). Extrasolar planetary orbits differ from the solar system’s in two ways: 
first is their semi-major axis distribution. Planets can orbit their parent stars much closer 
than the tenth of Mercury’s distance to the Sun. Second, at a median value of e = 0.27, their 
orbital eccentricities are so high that the only solar system analogs are small bodies such 
as asteroids and Kuiper belt objects whose orbits have been strongly stirred up by the gas 
giants. These differences are the more puzzling as it has been recognized that the distribution 
of semi-major axes and eccentricities of extrasolar planets resemble those of spectroscopic 
binary star systems, a finding that may hint to common formation or excitation mechanisms 
(Stepinski & Black 2001). 

Current theories of planetary formation suggest that Jupiter-like planets form in a disk 
of gas and dust orbiting a parent star. Producing planets with small semi-major axes is 
a natural outcome of the transfer of angular momentum between the disk and the planet 
(Ward & Hahn 2000). The overall effect of this transfer causes the planet’s orbit to shrink 
around the star. Radial migration can be stopped by the tidal interaction between the planet 
and the star (Lin et al. 2000) or through photo-evaporation of the disk’s portion exterior to 
the planet’s orbit and that is responsible for its inward motion (Hollenbach et al. 2000). 

The presence of large eccentricities is perhaps the most surprising feature of extrasolar 
systems as its unexplained origin contrasts with the known ability of the gas disk to force mi¬ 
gration. The transfer of angular momentum between the disk and the planet at the locations 
of their orbital resonances damps the planetary eccentricity so that the natural outcome of 
this interaction is a planet on a circular orbit (Goldreich & Tremaine 1981). This lead to the 
investigation of other possibilities: secular perturbations by distant companions (Holman et 
al. 1997; Terquem & Papaloizou 2002), resonant interactions within a multiplanet system 
(Chiang et al. 2002; Lee & Peale 2002), and encounters of a spatially extended multiplanet 
system with passing stars (Zakamska & Tremaine 2004). The advantages and drawbacks of 
these models have been reviewed by Tremaine & Zakamska (2004). The main recurring obser¬ 
vation is that these models provide a solution for some specific planetary system. Moreover, 
they do not explain the apparent uniqueness of the small eccentricities of the solar system’s 
planets, the occurrence of the largest eccentricities in some multiplanet systems for the outer 
more massive companions, and the origin of the similarity with spectroscopic binaries. In 
this regard, it is conceivable to assume that planets formed from a collapsing proto-stellar 
cloud but this would solve the problem only for the more massive known planets (Papaloizou 
& Terquem 2001). 

In this work, we aim to explain these observations within a single framework. We choose 
the premise of standard planetary formation theory that the natural outcome of planetesimal 
accumulation and gas accretion are planets on co-planar circular orbits. To account for the 
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similarity of the eccentricity distributions of extrasolar planets and spectroscopic binaries, 
we consider the excitation caused by physical processes that are weakly dependent on the 
local dynamics of the companion and which can therefore be modeled by an acceleration that 
does not depend explicitly on the relative position and velocity of the companion. In section 
2, we study the general properties of such an acceleration and show that stellar mass loss 
phenomena such as jets and winds are natural candidates for the excitation process. However, 
we do not specialize in the case of stellar jets because we prefer to present this theory as a 
phenomenological one. The reason is twofold: first, the current knowledge of stellar jets and 
winds both theoretical and observational is not sufficient to issue a definitive verdict on their 
role in the eccentricity excitation. Second, the dynamical mechanisms that arise because 
of such accelerations and that lead to eccentricity excitation and apsidal alignment are not 
concerned directly with the details of jet and wind generation. In section 3, we demonstrate 
the basic eccentricity excitation mechanism for the simplest finite duration acceleration, one 
whose direction is constant with respect to the star-companion system. In the remainder of 
the paper, we consider the various constraints on the basic acceleration model of section 3 set 
by: the acceleration’s precession (section 4), the radial migration of the companion (section 
5), the companion’s orbital stability and the consistency with the host star’s galactic motion 
(section 6), the secular perturbations from mutual gravitational interactions in multiplanet 
systems (section 7), and the presence of distant companions (section 8). The application of 
this theory to the solar system is further discussed in section 8. Section 9 contains concluding 
remarks. 


2. Phenomenology 

In this section, we discuss the form of the eccentricity generating processes that would 
apply equally to planetary and stellar companions. To this end, we seek processes that are 
weakly dependent on the local orbital dynamics of the companions. By this we mean that 
the interaction does not depend explicitly on the relative position and relative velocity of 
the star-companion system. With this in mind, we note that the simplest way to excite 
the eccentricity of a companion’s orbit is to subject it to a constant acceleration. Stable 
excitation of the eccentricity is possible if the constant acceleration is smaller than or com¬ 
parable to the gravitational accelerations between the star and its companions. Writing the 
perturbing acceleration as A = 4u where A represents its magnitude and u is a unit vector, 
the eccentricity excitation time is proportional to v/A, where v is the keplerian velocity of 
the companion. This determines two regions around the host star whether the ratio of the 
excitation time to the orbital period is larger than or comparable to unity. The former is 
located closer to the star as the excitation time decreases with the relative distance of the 
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star to the companion. We call this region the secular excitation region as the eccentric¬ 
ity increase builds up slowly after each rotation of the companion around the star. The 
eccentricity excitation in this case is derived in sections 3 and 4. The outer boundary of 
the secular excitation region is the sudden excitation region where the excitation time is 
comparable to or smaller than the orbital period. The dynamics in the sudden excitation 
region is considered in section 6. With the acceleration magnitudes that we determine below, 
the observed planetary companions fall into the secular region. Here and if the perturbing 
acceleration were strictly constant, the corresponding force would be conservative and the 
eccentricity excitation would simply result in periodic oscillations between the original value, 
zero, and some maximum value determined by the parameters of each individual system. We 
would require that A be finite for a duration comparable to or smaller than the oscillation 
period induced by the acceleration in order to ensure a finite residual eccentricity. In the 
secular region, the duration is naturally larger than the orbital period of the planet. In the 
sudden excitation region, the duration of excitation can be comparable to or smaller than 
the orbital period of distant companions and the outer minor bodies of the planetary system 
in order to ensure stability (see sections 6 and 8). We can obtain a minimal estimate of the 
acceleration’s maximal amplitude, A 0 , in the secular region by requiring that the eccentricity 
excitation time, be shorter than the age of the planetary system ~ 10 9 years. This yields: 
A 0 > 3 x 10~ 16 (u/10kms -1 ) kms“ 2 . 

The direction of acceleration, u, needs to be specified in an inertial frame as the accel¬ 
eration does not excite the eccentricity of a companion located in the secular region if it is 
perpendicular to its orbital plane. The knowledge of the current state of the solar system 
suggests two options: a local one related to the mean orbital plane of the planetary system, 
and a global one related to the motion of the system in the Galaxy. For the local option, 
the vector u is referred to the total angular momentum of the planetary system or the host 
star’s rotation axis. In the solar system, we would expect u to be closest to the angular mo¬ 
mentum direction because of the small eccentricities of Jupiter and Saturn. The direction of 
u can be made independent of the planetary planes if it is referred to the Sun’s rotation axis 
which is inclined by 6° to Jupiter’s orbital plane normal. On the scale of a planetary system, 
the simplest way to achieve an acceleration that is weakly dependent on the dynamics of 
the companion is by stellar and disk mass loss. The global option can be motivated by the 
galactic configuration of the solar system: the planets’ mean orbital plane is inclined to the 
Sun’s orbital plane in the Galaxy by the large value of 60 degrees. Jupiter and Saturn have 
small eccentricities compared to most of the known giant planets that are not located in their 
parent stars’ tidal zones. The direction of the acceleration can then be tied to that of the 
star’s motion in the Galaxy with the expectation that the largest eccentricities correspond to 
the situation where the planetary orbital plane is coincident with the parent star’s galactic 
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motion plane. The simplest two options for u are a constant vector that lies in the galactic 
mid-plane and does not depend on the star’s motion, and a vector directly related to the 
star’s motion such as its velocity in the Galaxy, v s , or the orthogonal thereof bxv 9 where 
b is a constant vector. On the larger scale of the solar neighborhood, there are no obvious 
processes that do not depend explicitly on the relative position and relative velocity of the 
star-companion system and that would excite the companion’s eccentricity. 

Jets and star-disk winds are ubiquitous features of star formation and star-disk inter¬ 
action (Eisloffel et al. 2000; Hartigan et al. 2000) and are natural candidates for the origin 
of the perturbing acceleration. Inferred mass loss rates for known young T Tauri stars lie in 
the range ~ 1O -8 M 0 year -1 to 1O -1O M 0 year -1 and may be two orders of magnitude larger 
depending on the way the rate is measured from the luminosity of forbidden lines (Hartigan 
et al. 1995; Hollenbach 1985; Kwan & Tademaru 1995). The mass loss process needs to be 
asymmetric with respect to the companion’s orbital plane in order to produce a residual 
acceleration. Interestingly, a number of bipolar jets from young stars (Hirth et al. 1994; 
Lavalley et al. 1997) are known to be asymmetric as the velocities of the jet and counter-jet 
differ by about a factor of 2. Mass loss processes in young stars therefore yield maximal 
accelerations: 


H 0 ~ 10 -13 


M \ 
1O -8 M 0 year -1 J 


300 km s 




( 1 ) 


where M is the stellar mass and v e is the outflow’s high velocity component. The implicit 
proportionality constant depends on the relative mass loss ratios M and ejection velocities 
of the jet and counter-jet. 

Integrated over the duration of acceleration, the star acquires a residual velocity that 
must be smaller than its galactic velocity else the star is ejected. The residual velocity is 
further bound by the known random component of the stellar motion. In section 8, we 
further discuss the mass loss rate values after we consider the constraints obtained from the 
companion’s orbital stability and from the eccentricity excitation of multiplanet systems. 
Another requirement for stellar and disk mass loss to be efficient in the secular region is that 
the momentum communicated to the star be inclined with respect to the companion’s orbital 
plane as it is the component of the acceleration that lies in the companion’s orbital plane that 
excites the eccentricity. Studies of molecular outflows suggest that some jets precess over 
timescales from 10 2 to 10 4 years (Davis et al. 1997; Eisloffel et al. 1996; Terquem et al. 1999). 
As most T Tauri stars are known to be in multiple systems, a possible way to achieve a jet 
precession is by warping the accreting disk’s plane through the gravitational perturbation of 
a stellar companion on an inclined orbit (Terquem et al. 1999). This would not imply that 
jet acceleration may only explain the large eccentricities of planetary companions in binary 
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systems. We show in sections 6 and 8 that the acceleration may force the outward migration 
and possibly the ejection of the stellar companion. 

In this paper, we do not specialize in a specific acceleration generating process and there¬ 
fore we do not enter into such specific details as asymmetric jet generation and precession 
by an accreting circumstellar disk tidally interacting with a stellar binary component. We 
study the eccentricity excitation by representing the perturbing process by a finite duration 
acceleration that is independent of the relative position and velocity of the system. In this 
respect, the application of this theory to the efficiency of stellar jets in exciting eccentricities 
is phenomenological and yields useful bounds on the star-disk mass loss rate. We point out 
that the excitation process that is modeled by this acceleration may depend implicitly on the 
relative position and velocity of the companion. For instance, the accretion that powers the 
stellar jet is driven by viscosity which can itself be driven by the planetary embryos through 
the shocking of the pressure waves they launch in the gas disk (Goodman & Rafikov 2001). 
In this case, the acceleration from the jet will depend implicitly on the companions’ local 
dynamics. 

The equation of motion is that of the two-body problem modified as follows: 


dv 

dt 


G(m + M) 


x + A(t) 


( 2 ) 


where A is the relative acceleration acting on the star-companion system, x and v are the 
relative position and velocity, M and m are the masses of the star and the companion. The 
acceleration is finite over a typical duration that we call r. For simplicity, we choose the 
amplitude of acceleration to be one of Ai{t) = A 0 H(t) exp —t/r and A 2 (t) = A 0 / cosh(t — 
to )/t where H(t) is the Heaviside unit step function and A 0 is the maximum acceleration 
amplitude. We will use the notation Aq when we refer to a constant amplitude acceleration. 
The connection between the perturbing acceleration and the star’s random motion in the 
Galaxy makes the residual velocity, V, a relevant parameter in the problem; it is defined as: 


V 



A(f) dt. 


( 3 ) 


For a constant u and the two specific forms we use, Vi = A 0 t and V 2 = nA 0 r. We note 
that it is likely that the companion is subject to radial migration due for instance to the 
tidal interaction with the circumstellar disk. The effect of migration on the eccentricity 
excitation is discussed in section 5 where we modify (2) to account for migration and show 
that our conclusions about the eccentricity remain valid in the presence of migration. By 
using equation (2), we also neglect the tidal effect of the Galactic potential; this is justified 
by the typical sizes of the planetary and stellar systems considered here. 
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The eccentricity excitation time, v/A a , allows us to split the excitation problem into two 
cases: the first is when the variation of the direction of A is small on the scale v/Aq in which 
case the direction of A can be held constant to a good approximation. The eccentricity 
excitation is demonstrated in this case in section 3. The second case is when the direction 
of A varies on a timescale comparable or shorter than v/Aq. This is the general case of 
precessing jets which is examined in section 4. 


3. Eccentricity excitation 

Under the effect of a constant-direction acceleration, A = A(t)u, the variation of specific 
angular momentum, h = x x v, is h = x x A. As u is constant, only the projection of the 
angular momentum along the direction of acceleration, h • u, is conserved. The rate of 
change of the angular momentum in the direction orthogonal to h and u is obtained from 
(uxh) h = A (u-x) (h-u) and shows that if u lies in the orbital plane, the orbit does not gain 
inclination with respect to its initial state. To find the eccentricity and inclination excitation 
of the companion’s orbit in the secular excitation region, we average the equations of motion 
and in particular the perturbing acceleration over the fast orbital motion of the companion. 
This is valid because the excitation time and the acceleration’s duration are larger than 
the orbital period of the companion. We can therefore assume that the acceleration is 
independent of time in deriving the secular equations. A constant acceleration derives from 
the potential R — A • x. Its average with respect to the orbital motion is derived in the 
Appendix and reads: 

3 3 

(R) = --aA(t) ■ e = --A(t) aesin(zu — Q) sin / (4) 

where e = v x h/G(m + M) — x/|x| is the eccentricity vector of magnitude e and a is the 
orbital semi-major axis. The last expression is obtained by choosing the 2 -axis along u; 
in this case, w, and / are the longitude of pericenter, the longitude of ascending node 
and the inclination of h with respect to u. The freedom of choosing a vector base in this 
plane results from the conservation of h • u and leads to the appearance of the combination 
u) = vd — fl, the argument of pericenter. 

The conservation of u - h can be written as a/ 1 — e 2 cos / = cos Jo, where Jo is the initial 
inclination of the circular orbit, and allows us to reduce the problem to a single degree of 
freedom with the potential: 



1 — e 2 


e sin a;. 


( 5 ) 



The eccentricity and pericenter evolution is obtained from 


• = _^LHZ • = (6 ) 

na 2 e du ’ U na 2 e de 

where n = y/G{M + m)/a 3 is the companion’s mean motion (Roy 1991). In a conservative 
system where A(t) = 4 is constant, e and u follow curves of constant (R). There are 
equilibria at c o — ±90° and e = \/2sin(/ 0 /2) corresponding to I = cos _ 1 (a/cos/ 0 ). The 
maximum value of e is sin I 0 and corresponds to the cycle of initially circular orbits. In 
Figure (1), we show the orbits of the conservative case for an inclination Iq = 30°. 

For initially circular orbits, the identity (R) = 0 leads to u = 0 or 180° throughout the 
evolution. The eccentricity equation then reads: 


3 A(t)e 
2 na 



(7) 


where e is the sign of coscn which is set by the requirement that e > 0. For initially circular 
orbits, the discontinuous changes of uj between 0 and 180° correspond to the crossing of the 
plane 1 = 0 during the oscillation of / between —I 0 and Jo; it is the result of the geometric 
requirement that 0 < I < 180° in standard keplerian variables. Finally, the eccentricity that 
is excited by the perturbing acceleration is: 


e(T) 



"3 r T . ' 


sin 

- / Ait) dt 

.2 na .4^ 

sin Iq 


( 8 ) 


The inclination is obtained from cos / = cos/o/a/1 — e(T) 2 . In the conservative case, A(t) = 
A 0 , e would oscillate between 0 and sin/ 0 at the excitation frequency: 


3 |4) 
na 


(9) 


Examples of such oscillations that were obtained from the direct integration of the full 
equations of motion are shown in Figure (2). The agreement between the secular solution 
and the results numerical integration is perfect and is due to the fact that A is independent 
of the relative position and velocity. To optimize the excitation of a finite eccentricity from 
an initially circular state, the duration of acceleration needs to be smaller than half the 
oscillation period: r < 7ma/3|4|- A less conservative criterion for eccentricity excitation 
is that the adiabatic condition, r 3> ima/3\A 0 \, that would ensure a long term decrease 
of e to zero is not satisfied. This allows the eccentricity to experience a few oscillations 
before settling down to a finite value. For multiplanet systems, this situation leads to strong 
interactions among the companions which can lead to ejections offering a possible explanation 
for the existence of many extrasolar planetary systems with one large planet. Examples 
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of eccentricity excitation at three different semi-major axes (i.e. three different excitation 
frequencies) are shown in Figure 3 for the two forms of the acceleration’s amplitude where r 
has been chosen to yield the same residual velocity V and hence the same final eccentricities 
for each semi-major axis. Figure 3 also illustrates the dependence of the excitation amplitude 
on the ratio of the duration to the excitation time. We note that because eccentricity 
excitation in the secular region is a slow process compared to the orbital time, the convolution 
of the dynamics under the conservative acceleration A 0 with a finite time window has the 
effect of shutting off the excitation at some eccentricity value depending on the duration. 

We conclude that the eccentricity e can be excited np to sin Jo and is largest if the initial 
orbital plane contains the direction of acceleration (Jo = 90°). As (R) = 0 for initially circular 
orbits, the argument of pericenter, cn, and the longitude of pericenter, hi, remain at zero. We 
show in sections 7 and 8 that this forcing of the pericentre to be aligned with the direction 
of acceleration favors apsidal alignment in multiplanet system. The inclination / decreases 
from its initial value pushing the orbital plane away from the direction of acceleration. We 
also note a interesting feature of this model that the final eccentricity increases outward so 
that in multiplanet systems, the farthest planets may have the largest eccentricities. This is 
a consequence of the relative strength of the gravitational acceleration and the perturbing 
acceleration as a function of distance. 


4. Precessing accelerations 

The inclination of the direction of acceleration u is crucial to the excitation of eccen¬ 
tricity in the secular region. If the acceleration is due to a stellar jet that precesses, we 
expect the excitation amplitude to depend on the ratio of the excitation frequency tia to 
the precession frequency f 1a and that resonant forcing is possible when a match occurs, a 
situation that is likely since tia is an increasing function of the semi-major axis a. In the 
following, we solve the conservative excitation problem with precession, that is we choose 
A = A 0 u(t) where u rotates at the rate Qa- For this problem, it is more convenient to 
derive the eccentricity vector and angular momentum evolution with secular perturbation 
theory in vector form (Milankovitch 1939) as: 

4 = -^k x V e (i?) +V k (i?), 

dr na- rur 

H1 1 

— = — kx V k (i2) + — ex V e (i?>, 

dr na na z 


( 10 ) 

( 11 ) 
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where k = h/ a JG(M + m)a is the non-dimensional angular momentum vector. Using the 
secular potential (5), these equations reduce to: 

de 32L 0 , dk 3A 0 

— =-k x u, — =- 

at 2 na at 2 na 

du 

dt 


e x u 


= n* x u, 


( 12 ) 

(13) 


where 12^ is the rotation vector associated with the precessing acceleration. With precession, 
the constants of motion are given as: 


3A, 


u • e + ft A ■ k = Ci, and 


3A n 


u • k + Q A ■ e = C r 2- 


(14) 


2 na 2 na 

To find the excitation under the precessing acceleration, we introduce the two vectors k-j_ = 
k ± e which decouple the equations of motion and yield: 

dk ± 32L 0 


uxk±, and k± ■ 


3A 0 

il A =F- u 

2 na 


= C .* 


(15) 


dt 2 na 

where C± — Cj ± C 2 . Note that |k±| = 1 is conserved, a consequence of k ■ e = 0 and 
e 2 + k 2 = 1. We seek a solution for the vectors k± through their projection on the basis 
made up of fl A , u, and ft A x u. Simple algebra shows that ft a • k± satisfy the equations: 

d\n A ■ k ± ) 


dt 2 


\p*a /4 T -F ^a(S^a ■ u)J {fl A ■ k±) — 

[(f A • u)/2] C±, 


(16) 


which are those of two harmonic oscillators of frequencies n±: 

n 2 ± = n 2 A /4: + Q a T n A (Q A • u) = {ft A t &k) 2 - (17) 

where ftk = 3Aou/2na is the instantaneous rotation vector of k + . The solutions of the 
previous equations are: 

[ 212(4 =F n A (fl A ■ u)] C± 


f2 4 • k+ = 


2n 2 


+ K± cos (n±t + <f>±) 


(18) 


± 


where K± and (j)± are constants to be determined from the initial conditions. The projections 
on u are obtained from the constants of motion as: 


=Fu • k ± = 


[n A ^2{Cl A -u)\C, 


± 


2A '4 


2 n± n A 


cos (n± t + <f>±) 


(19) 


The projections along fl A x u are found by noting that (f l A x u) • k± = ±2d(f2,4 
k±)/d t/n A which lead to: 

(Oi x u) • k± = =f— ^ ± sin (n± t + <f>±). 
n A 


( 20 ) 



For initially circular orbits that interest us, the integration constants are given by: 


K± cos (f>± 


K± sin <j)± 


n\ + 2n A {fl A • u)] C x 
4 n± 

[2n A {ft A • u) + 4 fl\] C 2 
4 ri± 


+ 


uaC 0 
2 n± 


( 21 ) 

( 22 ) 


where C' 0 denotes the initial value of (f2^ x u) • k±. The eccentricity and inclination expres¬ 
sions are obtained from e = \J{1 — k + • k_)/2 and cos/ = a 4 -(k+ + k_)/2fl A Vl - e 2 (note 
that the components of k± are not given in an orthogonal basis; to recover k + • k_ easily, an 
additional step is needed and consists of writing u in terms of its invariant part along fl A and 
its precessing component along the unit vector as u = V 1 - • u) 2 U 4 + (f2 A -u) ft A 2 fl A 

and substituting it into equations (19) and (20)). This completes the solution of the excita¬ 
tion by a precessing acceleration. 

The eccentricity excitation differs from that without precession in three ways: (i) the 
excitation amplitude depends on the ratio n A /2fl A , (ii) the motion involves two fundamental 
frequencies n± if ft A and u are not orthogonal, and (iii) large eccentricity excitation becomes 
accessible at all initial relative inclinations through the resonance n A = 2fl A where the 
frequencies n± in the denominators of the amplitudes of k± become small. 


In the following, we apply these findings to the case of a precessing jet whose rotation 
vector ft A is parallel to the initial angular momentum vector or equivalently the initial vector 
k. The constants of integration are given as: C 0 — 0, C\ — and C 2 = —n A cos(a )/2 where 
a is the angle between u and ft A . Denoting by p the frequency ratio n A /2ft A and using the 
solution derived above, we find: 


e = 


cos / = 


p 2 sin 2 a 
4 v\v 2 _ 


[2(3 + p 2 ) — 4(1 +pcosa) cosu + t — 4(1 — pcosa) cos v_t 


+ (1 — p 2 + i/ + i/_) cos(z/ + — v_)t +(1 — p 2 — v + v_) cos(z/ + + v_)t\ (23) 

~ 2 9 1 / -- ^ [( P 4 - P 2 + 2 + P 2 \p 2 - 3] cos 2a) 

2v+vi vl — e 1 

+p 2 sin a 2 (p 2 + 1 + 2 p cos a) cos u + t + p 2 sin a 2 (p 2 + 1 — 2 p cos a) cos z/_f](24) 


where v± = (n±/2fljC) 2 = p 2 + 1 + 2pcosa and t is normalized by 12 ^; for dehniteness we 
take Aq > 0 in what follows. 


The location where the frequency match, p — 1, occurs defines the nominal reso¬ 
nant semi-major axis a res = G(M + rn)(2flA/3Ao) 2 . The frequency ratio can be written 
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as p — \Ja/(ires- We consider the differences introduced by precession in the three re¬ 
gions: far inside resonance (p <C 1), far outside resonance (p 3> 1) and the resonance 
region p — 1. Far inside resonance (p -C 1), the jet precesses faster than the eccentric¬ 
ity excitation leading to a reduction of the eccentricity amplitude from since to 2p since 
as e = p\ since) [(3 + cos[2pcos(ce)t] — 4cos[t] cos[pcos(ce)t])/2] 1 ^ 2 . Far outside resonance 
(p> 1), the jet’s precession is slow compared to the eccentricity excitation so that the latter 
is described by the expressions given in section 3 for a constant-direction acceleration; the in¬ 
clination, however, is not. The slow precession causes a modulation of the inclination oscilla¬ 
tion between 0 and sin 2a as cos/ = [sin 2 a cos (pt) cos(cos[a]2)+cos 2 a]/[l— sin 2 a sin 2 (pt)]^ 2 . 
In the resonance region, the proximity of p to unity increases the denominators of the eccen¬ 
tricity expression (23) which leads to eccentricities close to unity. The width of the region 
around resonance where large eccentricity values are reached increases with the jet angle a. 
These features are illustrated in Figures (4) and (5) where we plot the expressions (23) and 
(24) for two jet angles a = 1° and 30°, an excitation period 2 tt/ua = 10 4 years, and the four 
values of p: 0.05, 0.5, 0.9, 1 and 50. Figure (4) also includes the result of the integration of 
the full equation of motion (2) and shows that the solution equations are indistinguishable 
from the unaveraged numerical solution. This agreement results from the fact that the ac¬ 
celeration is independent of the relative motion and the relative velocity. The figures show 
that the resonance region does not extend far around p = 1 for a jet angle a = 1°, as the 
maximal amplitude for p = 0.9 is 0.18 while for a = 30° the resonance region is much wider. 
We note that consistency in the use of the secular potential requires that p is not too small or 
equivalently that 12 a <C n. Numerical integrations of the full equations of motion show that 
the eccentricity and inclination expressions can be used as long as 12 a < 0.1 n. Moreover, as 
the eccentricity excitation time is tia, no resonant forcing occurs when 12 a = n in the secular 
region (tia -C n). 

Precession-driven resonant excitation therefore provides a possible way to raise the ec¬ 
centricity even in low jet angle systems. If a companion happens to be at or cross the 
excitation region, because of disk-driven migration, not only that the eccentricity will grow 
but the planet leaves the disk as its inclination is excited in phase with the eccentricity. The 
tidal interaction of the disk with the planet will not prevent the eccentricity excitation if the 
precession period of the jet, 10 2 to 10 4 years (Davis et al. 1997; Eisloffel et al. 1996; Terquem 
et ah 1999), is shorter than the viscous time of the disk, 10 5 to 10 6 years (Adams & Lin 
1993). Such a situation offers a possible prospect for stopping migration while exciting the 
eccentricity provided that (i) the disk is stable to the perturbing acceleration and that (ii) 
the planet’s migration time is longer than the excitation time. These results are discussed 
further after we constrain the magnitude of the perturbing acceleration in section 8. 



5. Effect of radial migration on eccentricity excitation 


The changes in the orbital semi-major axis a of the companion have so far been neglected 
because in the secular excitation region where tia is smaller than the mean motion n, a 
constant acceleration will produce small periodic oscillations of the semi-major axis with 
frequency n. It is possible that radial migration occurs indirectly even for an acceleration 
that is weakly dependent on the local orbital dynamics. Consider for instance the case 
where acceleration is caused by mass outflow from the star-disk system. The planet’s orbital 
revolution is determined by the matter content inside its orbit. If the total mass of that 
content varies over timescales larger than the orbital period the conservation of angular 
momentum leads to the radial migration of the companion but does not affect its eccentricity 
(Jeans 1924) leaving initially circular orbits invariant. This conservation leads to the relation 
Ma = Modo which in the case of mass loss induces outward migration. To estimate the 
related migration, we note that the disk’s mass inside the planetary companion’s orbit is 
at most several percent of M 0 . If all of this mass is ejected from around the star, the 
companion’s semi-major axis expands by a corresponding several percent showing that in 
this case migration is not significant. 

Disk-companion interactions usually yields an inward radial migration. To find the effect 
of migration on the eccentricity excitation, we simulate the disk-companion interaction by the 
addition of a Stokes-type drag — kw to the equations of motion (2) where k may be a function 
of time whose characteristic timescale is larger than the orbital period. Simple algebra 
shows that in the absence of external acceleration, the drag term conserves the modified 
angular momentum // h where rj = exp f kdt. In terms of osculating orbital elements, the 
previous relation implies that a(l — e 2 ) r/ 2 is conserved showing that the drag term leads 
to orbital migration of the companion with respect to the star. The time dependence of k 
can be interpreted as the mathematical representation of the decay law of a; for instance, if 
a = ao{t/r a + l) _a , then k = —a(t/r a + l) _1 /2r a with nr a -C 1 where r a is the migration 
timescale that depends on the parameters that regulate angular momentum transfer between 
the companion and the gas disk. A constant k implies an exponential decay of a. When the 
migration time is larger than the companion’s orbital period, which is usually the case, the 
average effect of this drag term conserves the orbit’s eccentricity and the planet stays in its 
initial circular orbit while migrating with respect to the star. 

Radial migration would therefore affect the eccentricity excitation only through the vari¬ 
ation of the excitation frequency tia which is where the semi-major axis enters the excitation 
mechanism. This is illustrated in Figure (6) where we show an example of the excitation of 
eccentricity during migration. For multiplanet systems, the eccentricity excitation can be 
affected by the migration significantly because the acceleration has to compete with mutual 
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planetary perturbations that tend to precess the planetary orbits at rates that may be faster 
than the acceleration’s excitation frequency. 


6. Keplerian boundary and sudden excitation 


In order to be stable, a companion’s orbit needs to receive greater acceleration from the 
star than from the perturbation. This condition delineates the keplerian region around a 
star as that for which G(M + m)/r 2 > |A|. The natural limit for small orbital perturbations 
is located closer to the star where the frequency n A becomes comparable to the local mean 
motion n of the companion. Near this limit, the forced periodic oscillations of the semi-major 
axis a are reinforced by the eccentricity and acquire large amplitudes. Denoting by ak p i r the 
semi-major axis of the keplerian boundary where tia = n, we obtain an expression for the 
magnitude of the perturbing acceleration as: 


|d 0 | ~ 2 x 1(T 12 
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which for a^pir = 10 3 AU is an order of magnitude larger than the acceleration of the solar 
system in the Galaxy, the match occurs near a^ph- = 3000 AU. The excitation period of this 
acceleration, Ta = 2n/n A , is given as: 
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When the boundary ak p i r of a given acceleration is located beyond 10 4 AU, the Galactic 
potential’s tide becomes important (Heislcr and Tremaine 1986) and must be included to 
determine the limits of the boundary. In this work, we will be concerned with smaller values 
of akpir- To illustrate the motion near the keplerian boundary, we show an example of an 
escape orbit of a conservative constant-direction acceleration with ak p i r = 10 2 AU and an 
inclination I 0 = 30° (Figures 7 and 8 ). The orbit’s initial semi-major axis is 68.5 AU. The 
characteristics of motion are not strictly keplerian as the companion hovers above the star. 
The semi-major axis shows periodic oscillations around 100 AU with a significant amplitude. 
These characteristics depend to a certain extent on the direction of acceleration. A study 
of the types of motion near the boundary is interesting but not central to the problem 
of the eccentricity excitation. We remark that the escape orbits offer an interesting way 
to expel planets from around their parent stars or equivalently to disrupt a binary stellar 
system. If a companion is formed near the keplerian boundary or is pushed out to it by 
a possibly remaining inner disk that followed photo-evaporation (Hollenbach et al. 2000; 
Veras & Armitage 2004), it could become unbound. Systems where the acceleration keeps 
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an approximately constant direction during eccentricity excitation loose companions without 
a risk of catastrophic encounters with the star if the initial direction of acceleration does not 
lie in the companion’s orbital plane such as that of Figure ( 8 ). 

We remark that the value of dk p i r corresponds to the innermost location of the stability 
boundary because the acceleration’s decay decreases tia and pushes the stability boundary 
outward. The acceleration’s finite duration may extend the keplcrian boundary depending 
on the ratio of r to the excitation time at the keplerian boundary of the conservative problem 
Ta{o> k P ir)/2. When r > T^(ak p i r )/2, orbits beyond ak p i r have enough time to acquire sufficient 
momentum to escape the gravitational pull of the star. When r < Th(dk p i r )/2, the stability 
region extends beyond ak p i r and its boundary is given by the semi-major axis where r ~ 
Ta(o> 0 o)/2 which is larger than dk p i r since Ta is a decreasing function of the semi-major axis 
a. Using the expression of T 4 and the residual velocity V, a <*, ~ G(M + m)U~ 2 , the location 
where the keplerian velocity matches V. 

We now determine the features of orbital excitation in the keplerian boundary for an 
acceleration whose duration, r, is smaller than the excitation time at the keplerian boundary. 
We consider the outermost orbits in the sudden excitation region for which the excitation 
time is very small compared to the orbital time. To these orbits, the sudden excitation 
imparts a near instantaneous velocity V. The changes in the orbital elements are found by 
expressing three conservation relations. First is the conservation of the potential energy as 
the position x is left invariant during the excitation. This is written as: 
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where dj, df, V\ and Vf are the initial and final semi-major axes and velocities of the companion. 
Eliminating the potential energy term leads to: 
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For an initially circular companion orbit in a coordinate system where the z-axis is chosen 
along V (section 3), the previous equation becomes: 


2U sin/p cos 6 _ V 2 fi?Q) 

Of Oi yjG(M + m)di G(M + my 

where 6 is the longitude of the companion along its orbit, and Jo is the inclination of the 
orbital plane with respect to the direction of the residual velocity V. This equation deter¬ 
mines, for a given V, the final semi-major axis df as a function of d; and 6. We note that 
df can be larger or smaller than d; because of the inclination term in the energy equation 
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(29). This leads to an inward or outward migration of the companion. Note that when V is 
perpendicular to the initial orbital plane (Jo = 0). migration is always outward. 

The second relation is obtained from the conservation of the linear momentum, Vf = 
V; + V. For initially circular orbits, the projection of the velocity on the position vector 
yields: x ■ Vf = x • V. In the cases that we consider below where the companion is initially at 
the orbit’s nodes, x-Vf = 0 implying the conservation of the pericenter distance af(l—ef) = a x 
if migration is outward and the apocenter distance af(l + ef) = O; if migration is inward. 
This leads to the following relation that is valid at the nodes: 

ef = |1 — Oi/af|. (30) 


Elsewhere on the orbit, one needs to express the projections of the linear momentum conser¬ 
vation equation in order to derive the final eccentricity, longitude of pericenter and longitude 
of ascending node as functions of the longitude 9. The third relation comes from the conser¬ 
vation of the projection of the angular momentum on V: 


y Of(l — e 2 ) cos/f = y/ai cos J 0 , (31) 

where Jf is the final inclination of the companion. There are three particular initial semi- 
major axes that characterize migration: first is a esc , the semi-major axis beyond which 
the companion can escape at certain longitudes. Second is a^, the semi-major axis beyond 
which all companions are lost. Third is a out , the semi-major axis beyond which only outward 
migration occurs. Using equation (29), we find: 
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Migration is inward in the longitude range around 9 = 180° defined by: 
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Between a esc and a oo, escape occurs in the longitude range around 9 = 0 defined by: 


, qf 1 - V 2 /G(M + m) 
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(36) 


Figure (9) shows the final semi-major axes and eccentricities at two different inclinations 
I 0 = 0 and 20° for an acceleration corresponding to ak p i r = 300 AU and a duration r = 500 
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years resulting in a residual velocity V = 0.35 km s -1 . The companions were started at the 
descending node for J 0 = 20° to illustrate inward and outward migration. Also shown in 
the same figure are the results of the integration of the full equations of motion (2) with 
the form Ai(t) confirming the validity of the analytic expressions in the outer keplerian 
boundary. An example of migration with the form A 2 is used for the eccentricity excitation 
of the v Andromedae system in section 8.1. The dependence of migration on the orbital 
longitude is illustrated in Figure (10) for the inclined orbits with J 0 = 20°. The migration 
and excitation near the keplerian boundary offers a possible way to transport minor bodies 
in the solar system (section 8.2). For planetary companions, the radial migration that results 
from the interaction with the disk (section 5) has to be added to the migration in the sudden 
excitation region in order to ascertain the dynamics in the outer keplerian boundary. 


The previous study of orbital stability points out similar constraints for the host star 
that is accelerated in its motion within the Galaxy. Stellar jets for instance modify the 
star’s velocity but its galactic motion precludes that A 0 be larger than that the galactic 
acceleration for an arbitrarily extended time else the star is ejected. The acceleration’s 
duration, r, is therefore an important parameter that determines the star’s orbital stability. 
A limit on r can be set by requiring that the residual velocity, V, that contributes to the 
stellar random motion be smaller than the known velocity dispersion, ( v g ). The condition 
V < (v g ) yields: 


T < 10 5 3 x 10 12 km s 2 (v g ) 
A 0 10 km s~ 

which in terms of the size of the keplerian boundary is: 
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where we used V ~ A 0 r. For a precessing acceleration, the duration is extended by a factor 
1/ cos a: where a is the jet angle. In addition, the previous estimate is meaningful for a single 
acceleration event and does not account for an arbitrary time evolution of the angle a. 


7. Mutual planetary perturbations 

In section 4, we showed that the precession of the perturbing acceleration affects signif¬ 
icantly its ability to excite the companion’s eccentricity. This is also true if the acceleration 
does not precess but the system contains multiple companions. In this case, mutual gravi¬ 
tational interactions cause the eccentricity vectors to precess and possibly to be locked into 
resonance. We expect for instance that if the precession rates are much faster than the ex¬ 
citation frequencies, the maximum eccentricities will be reduced. To illustrate the effects of 
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mutual interactions, we model the Jupiter-Saturn system under the influence of a constant 
direction acceleration using a direct integration of the equations of motion. We assume that 
the planets formed in initially coplanar circular orbits and that the planetary orbital plane is 
inclined by 30° with respect to the direction of acceleration. The semi-major axes and masses 
are taken as aj = 5.20 AU, as = 9.55 AU, mj = 9.55 x 10 - 4 M o and ms = 2.86 x 1O - 4 M 0 . The 
unperturbed eigenfrequencies of the isolated Jupiter-Saturn system are (Murray & Dermott 
1999): 9.6 x 10 -4 ° year -1 and 6 x 10 -3 ° year -1 for e and —7x 10 -3 ° year -1 for I corresponding 
to the eigenperiods T) = 375 000 years, T 2 = 60 000 years and T 3 = 51 000 years respectively. 
Note that the eccentricity amplitudes have a period of 1/(1/60 000 — 1/375 000) ~ 71000 
years and that there is only one inclination eigenmode due to the freedom of choosing a refer¬ 
ence plane for Jupiter’s orbit. The evolution associated with mutual planetary perturbations 
with no external acceleration is depicted in the top panel of Figure (11) where we have as¬ 
cribed to the two planets their current eccentricities; the eigenfrequencies of the system are 
independent of e and / if these are small. We assess the effect of mutual perturbation by ap¬ 
plying two different acceleration strengths, Ao = 2 x 10 - 12 kms -2 and Aq = 2 x 10 - 14 kms -2 , 
corresponding to the two values ak p i r = 10 3 AU and ak p i r = 10 4 AU for the location of the 
keplerian boundary. Taking ak p i r = 10 3 AU yields T AJ = 4.4 x 10 5 years, T AS = 3.2 x 10 5 
years; a simulation of this configuration is shown in the bottom panel of Figure (11) for 
which we removed mutual interactions. The case where both mutual interactions and the 
perturbing acceleration are turned on is shown in Figure (12). The eccentricities of both 
planets are excited to 0.3 except that Saturn’s is slightly smaller than Jupiter’s. This can 
be explained by the fact that TJj ~ 7\ which implies a smaller eccentricity reduction for 
Jupiter. In the case of Saturn, we have T AS 3> T 2 which implies a stronger reduction than 
observed but the gravitational interaction with a more massive Jupiter compensates that ef¬ 
fect and forces a larger eccentricity. The dominant eccentricity oscillation period is modified 
to 2.6 x 10 5 years. For a kp i r = 10 4 AU, TJj = 4.4 x 10 6 years, T AS = 3.2 x 10 6 years. These 
timescales are much larger that the unperturbed eigenperiods implying a strong reduction 
of the excited eccentricities. Figure (13) shows that both eccentricities are smaller than 
0.01. The dominant eccentricity oscillation period in this case is modified to 3.6 x 10 5 years. 
We can find the minimum strength of the constant-direction acceleration that can excite 
Jupiter’s eccentricity to its current value 0.05. Numerical integrations yield ak p i r = 2600 AU 
corresponding to A 0 = 3 x 10 - 13 kms -2 . 

We remark that during the eccentricity and inclination excitation, apsidal alignment is 
maintained with the libration of fls — Uj and oos ~ oq around zero. This is the result of the 
forcing of the pericenters to be aligned with the direction of acceleration that we encountered 
in the basic two-body problem of section 3. Note however that the inclination oscillation 
about the direction of acceleration has been largely suppressed. For a,k p i r = 10 3 AU (Figure 
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12), the two planets’ inclinations oscillate around 26° with an amplitude of 3° and a mutual 
inclination of 0.2°. We can use this feature as another constraint on the initial strength 
of the constant direction force by imposing the current mutual inclination of Jupiter and 
Saturn of ~ 1°. Numerical integrations yield ak p i r = 590 AU and A 0 = 5.7 x 10 -12 kms -2 
with maximal eccentricities for Jupiter and Saturn of 0.45. 

We can obtain an additional constraint on the perturbing constant-direction acceler¬ 
ation by deriving its duration. We do this for the form A 1 with the following simplified 
configuration: ak p i r = 590 AU; the eccentricity excitation occurs during the first oscillation 
cycle; disk driven radial migration is neglected; and Jupiter and Saturn are assumed to have 
formed at their current locations. In this case, equation (8) leads to: 
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where e p is the current eccentricity for Jupiter and T e ~ 1.4 x 10 5 years is the excitation 
period obtained with mutual perturbations. For e p = 0.05, r ~ 5000 years and the residual 
velocity V ~ lkms -1 . The application of the model to the solar system is further discussed 
in section 8.2. We finally comment that some of the dynamical features that we have shown 
in this section can be reproduced analytically by applying the Laplace-Lagrange secular 
perturbation theory at small e and I, and combining its second order secular gravitational 
potential with the linearized acceleration potential (4) (i.e. substituting / for sin / in that 
expression). This would yield the linearized eigenfrequencies and eigenvectors of the system. 
We do not follow this approach because we are interested in large eccentricities and the 
possible effects of mean motion resonances and close encounters. 


8. Discussion 

The eccentricity excitation theory that we have presented is based on the fact that if 
a star-companion system is subjected to a relative acceleration that is weakly dependent 
on the local dynamics, large eccentricities can be achieved. The theory has the following 
parameters: 

1 . The magnitude of acceleration, A, which determines the excitation frequency of a 
star-companion system, and the extent of the keplerian region around the main star. 

2. The direction of acceleration, u, (unit vector) with respect to the initial orbital plane 
which determines the maximum eccentricity of a star-companion system as the projec¬ 
tion of u on the orbital plane. 



3. The rotation vector, ft a, which describes the acceleration’s precession whose effect on 
the eccentricity amplitude and resonant excitation is described in section 4. 

4. The duration of acceleration, r, which for a given eccentricity amplitude (larger than 
required) permits the selection of the final eccentricity value. 

When these parameters are specified, the final orbital state of a planetary or stellar com¬ 
panion can be selected. The orbital configurations of star systems with a single companion 
are easier to explain than those with multiple companions because the absence of additional 
perturbers relaxes the constraints on the acceleration’s magnitude or equivalently the size of 
its keplerian boundary Ok p i r - Since stellar jets and star-disk winds are ubiquitous features of 
star formation and since we have established that they can be dynamically responsible for 
the eccentricity excitation, the known similarities and differences of the eccentricity distribu¬ 
tions of planetary companions and spectroscopic binaries can be attributed to the similarities 
and differences of the physical environments that give rise to the accelerating mass loss pro¬ 
cesses and not directly to the companions’ formation or their dynamical interaction with 
the gas disk. In this regard, it would be interesting to be able to associate, for stars of 
a similar spectral type, similar residual velocities that reproduce the observed companion 
eccentricities. 

The excitation theory is built on a minimal assumption that the perturbing acceleration 
is explicitly independent of the local dynamics -the basic excitation mechanism (section 3) 
relies on a constant acceleration that is applied for a finite duration. As a result, multiplanet 
systems as well as planetary systems of binary stars put strong constraints on the accel¬ 
eration’s parameters. Once these parameters are specified, a single relative acceleration is 
applied to all companions, planetary and stellar alike, and sets the fundamental excitation 
frequencies and amplitudes of the system. These frequencies are influenced by the mutual 
perturbations of the companions as shown in section 7. In particular, the larger the compan¬ 
ions’ number, the more constrained the acceleration’s parameters. In section 8.1, we illustrate 
this situation with the v Andromedae binary system which contains multiple planets and 
a distant stellar companion whose observed location strongly constrains the acceleration’s 
strength. The application of this theory to the known sample of multiplanet systems is 
beyond the scope of this paper but it is ultimately the best way to ascertain whether the 
acceleration mechanism is responsible for the eccentricities of extrasolar planets. Two ef¬ 
fects need to receive particular attention: the companions’ possible radial migration which 
changes the fundamental excitation frequencies as seen in section 5, and the companion’s 
secular perturbations which circulate the eccentricities. Radial migration can be prescribed 
by the standard disk-planet interaction theories (Ward & Hahn 2000). The effect of secu¬ 
lar perturbations may be disentangled by determining for a given multiplanet system the 
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amplitudes of the secular modes from observed eccentricities and seeking the acceleration 
parameters that reproduce the amplitudes of the secular modes and not the observed eccen¬ 
tricities. This procedure would be useful for similar mass planets because the instantaneous 
eccentricities can be different from the individual secular mode amplitudes. The combination 
of both migration and secular perturbations relaxes some of the constraints on the model 
because of the additional freedom in the choice of the frequencies of the secular modes. Such 
an analysis may help elucidate the eccentricity-mass correlations of extra-solar planets. In 
the following, we further illustrate how to constrain the acceleration mechanism in the case 
of the v Andromedae system (section 8.1) and discuss the relevance of these results to the 
solar system (section 8.2). 


8.1. The v Andromedae binary system 

Three planets orbit v Andromedae (Butler et al. 1997), two of which have their apsides 
aligned (Rivera & Lissauer 2000; Lissauer & Rivera 2001; Chiang et al. 2001). The most 
recent observations of this system (Naef et al. 2004) yield the following orbital elements 
Ofe = 0.059 AU, eb = 0.020, c Ob = 241°, nib = 0.75 mj, a c = 0.821 AU, e c = 0.185, u c = 214°, 
m c = 2.25 mj, = 2.57 AU, = 0.269, u>d = 247°, and rn c i = 2.57 mj where the masses are 
line-of-sight values, the stellar mass is M = 1.3 M & . 

Using the numerical integration of the full equations of motion (2) to model the ec¬ 
centricity excitation with a constant-direction acceleration applied to the planetary orbits 
with their current semi-major axes but with initially circular co-planar orbits, we find that 
mutual planetary perturbations are strong enough to damp the excitation if the acceleration 
is smaller than A 0 ~ 10~ n kms -2 . The equivalent keplerian boundary is at ak p i r ~ 500 AU. 
Below this value, the present configuration can be recovered along with the apsidal alignment 
of the outer two planets. However, the stellar companion to v Andromedae (Lowrance et al. 
2002 ) puts an additional constraint on the excitation mechanism. At a projected distance of 
750 AU, the presence of this O.2M 0 companion leaves us two options: either the excitation 
by acceleration is ruled out because the companion lies far outside the keplerian boundary of 
the weakest possible acceleration that reproduces the planets’ eccentricities, or that it was 
initially inside the boundary and migrated by the sudden excitation that we discussed in 
section 6 . We point out that the projected distance of 750 AU does not translate necessarily 
into a semi-major axis as the companion’s orbit is likely to be eccentric. 

To test the second option, we choose an acceleration of the form A 2 with a conser¬ 
vative keplerian boundary at 300 AU corresponding to Aq ~ 3 x 10 _ 11 kms~ 2 , a duration 
r = 2000 years and an equivalent residual velocity V = 5.6kms _1 . The planetary angu- 
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lar momentum is inclined by 60° with respect to the direction of acceleration. The stellar 
companion’s initial orbit has a semi-major axis a = 298 AU, an eccentricity e = 0.3 and an 
inclination of 10 ° with respect to the direction of acceleration. At this large semi-major axis, 
the stellar companion’s perturbation of the planets is negligible; in particular, the eccen¬ 
tricity excitation by the Kozai mechanism (Holman et al. 1997) is not efficient because the 
corresponding excitation time (~ 10 ' years) is much larger the duration of acceleration and 
the eccentricity secular frequency of the isolated two-planet system (~ 7000 years). We have 
chosen such an eccentric and inclined stellar orbit to decouple the problem of eccentricity 
excitation for the planetary and stellar companions. We find that the acceleration produces a 
configuration similar to the observed one with stellar orbital elements: e = 0.5, a = 600 AU, 
/ = 50° (Figure 14). In particular, we remark that the apsidal alignment that is generic in 
the basic mechanism of section 3 for two bodies is maintained in the v Andromedae system. 
We also note that only when the acceleration’s strength is near maximum and the keplerian 
boundary nears 300 AU, does the orbit acquire a larger eccentricity. The numerical integra¬ 
tion also confirms the absence of the Kozai eccentricity perturbations in the orbits’ evolution 
as suggested by the timescale analysis. 

The acceleration’s magnitude of A 0 ~ 3 x 10 -11 kms -2 corresponds to a jet-driven mass 
loss rate M ~ 1O - 6 M 0 year -1 . This value is within the two orders of magnitude uncertainty 
in mass loss estimation (section 2). However, the inferred rate is a maximal value in time 
(e.g. Figure 14) and cannot be compared directly to the observed steady state values which 
are naturally smaller. An additional constraint on the excitation by jet acceleration is that 
the planet must not lie far inside the resonance location. In terms of the parameters we 
introduced, the resonant semi-major axis of the precessing jet is given as: 
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where T prec = 2 tt/VLa is the jet’s precession period. This estimate shows that the most effi¬ 
cient accelerating jets precess with periods of 10 4 years or larger. If the precession angle is 
large as it is the case here, then the large eccentricity region is extended inside the resonant 
semi-major axis by a sizable factor (e.g. Figure 5 with a = 30°). At a = a res /100, maximum 
amplitude is 0.17 for a = 60°. This simple example allows us to conclude that stellar-disk 
mass loss is a possible process to provide acceleration. Further modeling of jet accelera¬ 
tion with a more realistic amplitude time variation is needed to constrain the excitation 
mechanism. 
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8.2. The solar system 

The solar system was used in our model to motivate the possible options of the direction 
of acceleration using the small eccentricities of Jupiter and Saturn. The detailed analysis 
of the excitation by acceleration pointed out other ways to limit the eccentricity growth: 
to the maximum eccentricity set by the inclination of the acceleration with respect to the 
orbital plane, we add the amplitude reduction associated with a precessing acceleration and 
that associated with a weak acceleration (ak p i r > 2600 AU with Jo = 30° constant) that 
competes with mutual planetary perturbations. These options show that it is more difficult 
to identify the evolutionary track of low eccentricity systems. Modeling the migration of 
the solar system’s planets and its effects on minor bodies under a perturbing acceleration 
may offer a way to discriminate between the available options. The current state of the 
giant planets’ orbits is believed to have evolved from a much more compact system. Studies 
of the origin of the eccentricity and inclination distributions of Kuiper belt objects suggest 
that the four giant planets were confined between 5.5 AU and 13.5 AU at the time where the 
Kuiper belt was made up of a low eccentricity and inclination planetesimal disk (Malhotra 
1995; Gomes 2003). For such compact configurations, mean motion resonance crossing may 
produce eccentricities comparable to the current values for Jupiter and Saturn (Tsiganis et 
al. 2004). Regardless of the specifics of such mechanisms and if we admit that accelerating 
processes of the type studied here excited companion eccentricities for solar-type stars in the 
solar neighborhood, the same processes must have applied to the solar system and shaped 
it to a certain extent. In the context of acceleration by stellar jets and star-disk winds, the 
inclination of Jupiter’s angular momentum vector by 6° with respect to the Sun’s rotation 
axis is an indication that Jupiter’s orbital plane has had to evolve from an equatorial accretion 
disk or that the planet had formed in a warped disk. Both possibilities are favorable to the 
excitation by acceleration model. 

An interesting application of our model that does not require modifying the standard 
picture of how the solar system and its minor body populations have evolved is to choose 
the direction of acceleration nearly perpendicular to the planetary orbital plane. The accel¬ 
eration, however, may not be particularly weak thereby allowing smaller values of ak p i r - The 
main consequence of such an acceleration is to excite the eccentricities of bodies that lie near 
the keplerian boundary, eject those that lie far outside it, thus truncating the planetesimal 
disk, and outward transport those in its vicinity through the migration that we discussed 
in section 6 and illustrated in Figures (9) and (10). Migration could enhance the delivery 
of minor bodies to the Oort Cloud and explain the transport of Kuiper Belt outliers 2000 
CR105 and Sedna (90377) that, with perihelia larger then Neptune’s semi-major axis, still 
elude dynamical explanation (Brown et al. 2004; Morbidclli & Levison 2004). 



9. Concluding remarks 


In this paper, we have worked towards building a theory for the origin of the eccen¬ 
tricities of extrasolar planets. Proceeding with minimal assumptions, we showed that the 
planetary eccentricities can be caused by relative accelerations that depend weakly on the 
local dynamics. We have thus reduced the problem of the origin of eccentricities to the identi¬ 
fication of the physical processes that cause such accelerations. Possible processes are stellar 
jets and star-disk winds that accelerate the host star with respect to the companion and 
excite the random component of its galactic velocity. The origin of extrasolar eccentricities 
can therefore be related to the the random motion of disk stars in the Galaxy. 

The model has further applications to the dynamics of extrasolar planets and the solar 
system. If the duration of acceleration is larger than the excitation time at some semi-major 
axis, planets exterior to this radius will achieve maximum eccentricity and risk ejections by 
close encounters (section 3). This offers a possible explanation for single planet systems with 
large eccentricities. Rogue planets can also be produced if planets are pushed out to the 
keplerian boundary by disk-planet interaction (section 6). For the solar system, besides the 
possible eccentricity excitation, acceleration leads to the transport of minor bodies to the 
outer Kuiper Belt and the Oort Cloud. 

The excitation by acceleration model can be applied in other contexts beside the stellar 
and planetary and companions. The star-disk mass loss mechanisms, that can be responsible 
for the eccentricity of extrasolar planets, impart an equal acceleration to the protoplanetary 
disk. A precessing acceleration derives from the perturbing potential R = A-x = rA r cos (0 — 
6a) + zA z where 9a, A r , A z are the acceleration’s components in a cylindrical coordinate 
system referred to the disk. This potential could be able to excite the rn — 1 slow modes 
of the disk (Tremaine 2001). Such modes are interesting because of their large wavelength 
that can be comparable to the size of the system. For instance, the disk could develop rigid 
precession if the speed of sound waves is larger than the precession rate similarly to the case 
of tilted disks perturbed by a binary component (Papaloizou & Terquem 1995). The feedback 
of a jet-generated precessing acceleration on the disk may have important consequences for 
sustaining the acceleration. 

The excitation mechanism also applies to galactic dynamics. Large-scale wind phe¬ 
nomena (Bland-Hawthorn & Cohen 2003; Shopbell & Bland-Hawthorn 1998) accelerate the 
Galaxy and alter the galactic rotation curve by decreasing the circular velocity. The analysis 
of this paper can be modified to ascertain the efficiency of galactic winds in the excitation 
of the stellar random motion by substituting the galactic potential for the stellar potential 
in equation (2). Further applications of galactic winds include the stability of large scale 
structures in disk galaxies, the onset of m — 1 elliptic distortions, and the effect of the 
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change in the Sun’s Hill sphere in the Galaxy on the dynamics of the Oort Could. 
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A. Averaged potential of a constant acceleration 

In the following, we derive the averaged potential of a constant acceleration A. The 
acceleration (or equivalently its force) derives from the potential R = A • x where x is the 
position vector. The averaged potential is (R) = A - (x) and therefore only the average of the 
position vector x needs to be calculated. The companion’s unperturbed orbit is keplerian 
and in an orthonormal basis can be written as: 

x — r (cos hi cos(cn + /) — sin 72 sin(cn + /) cos I ), (Al) 

y — r (sin hi cos (a; + /) + cos 72 sin(a; + /) cos 7), (A2) 

z = r sin (a; + /) sin/, (A3) 

where the radius is r = a(l — e 2 )/(l + ecosf) and the orbital elements a, /, e, u,I, and 
72 are respectively, the semi-major axis, the true anomaly, the eccentricity, the argument of 
pericenter, the inclination and the longitude of the ascending node. The time average can 
be replaced with a true anomaly average by using the conservation of angular momentum 
d t/P = r 2 d//27ra 2 \/l — e 2 , where P is the period. Applying this integration to the position 
vector gives: 

(x) =-|ae-(/= 0). (A4) 

2 r 

The direction of x at / = 0 is that of the eccentricity vector (or Runge-Lenz vector), 
e = v x h/G(M + m) — x/r. The averaged potential is therefore: 

(R) = —- a A • e. (A5) 

Inspection of the position vector equations shows that the simplest expression corresponds 
to choosing the x-axis as the direction of acceleration. This leads to (R) = — 3aeAsin(zu — 
Q) sin 7/2 where w = u + is the longitude of the pericenter, the conjugate variable of the 
eccentricity, e. 
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Fig. 1.— Contour plots of the acceleration potential (5) in the eccentricity e and argument 
of pericenter cn(°) plane. The direction of acceleration makes an angle J 0 = 30° with respect 
to the companion’s angular momentum vector. The time evolution of the two orbits (e = 0, 
u = 0) and (e = 0.3, u = 90°) is shown in Figure (2). 
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Fig. 2.— Time evolution under a conservative acceleration. The eccentricity e, argument 
of pericenter cu(°) and inclination J(°) are shown for an initially circular orbit e = 0 (solid) 
and an orbit librating about the secular resonance lo = 90° with an initial eccentricity 
e = 0.3 (dashed). The semi-major axis is identical for both orbits and is set to unity. The 
acceleration corresponds to a period of 10 4 years at 1AU. The plots were obtained by the 
numerical integration of the full equations of motion (2). 
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Fig. 3.— Eccentricity excitation by time-dependent constant-direction accelerations. The 
full equations of motion (2) are integrated numerically with the forms A\ (two upper panels) 
and A 2 (two lower panels) and Iq = 30°. The parameters are: Ao = 2.21 x 10 -11 kms“ 2 
and to = 2.3 x 10 4 years. The oscillation period at 1AU is 1.11 x 10 5 years. The timescale 
r is chosen so that V =5kms _1 ; r = 7200 for Ai and r = 2300 years for A 2 . The curves 
correspond the semi-major axes: 1AU (solid), 32 AU (dashed) and 128 AU (dotted). Note 
how the final e and / are equal under the two different accelerations at each semi-major axis. 
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Fig. 4.— Excitation of the eccentricity (solid) and inclination (dashed) by a nearly perpen¬ 
dicular precessing jet with an angle a = 1°. The companion’s orbit is located at a = 1 
AU and evolves from a circular orbit in a plane orthogonal to the jet’s precession axis. The 
acceleration is Aq = 2 x 10 -lo kms~ 2 yielding an excitation time of 2 tt/tia = 10 4 years. From 
top to bottom, the panels were obtained from equations (23-24) with the frequency ratios, 
p: 0.05, 0.5, 0.9, 1 and 5 - the precession period is 2 p x 10 4 years. The symbols correspond 
to the numerical integration of the full equations of motion (2) and show that the agreement 
with the averaged analytical solution is perfect. 
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Fig. 5.— Same as Figure (4) for a jet angle a = 30°. From top to bottom, the panels 
correspond to the frequency ratios, p: 0.05, 0.5, 0.9, 1 and 50. 
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Fig. 6.— Effect of radial migration on eccentricity excitation. In this numerical integration of 
(2) with a constant direction acceleration of magnitude Ao = 2 x 10 -10 kms -2 and inclination 
Jo = 30°, a Stokes drag term, — kv with /c -1 = 36 500 years, was added to the equations of 
motion. From top to bottom, the panels show the eccentricity, argument of pericenter, the 
inclination and the semi-major axis. The companion is located at 5.2 AU on an initially 
circular orbit. The initial excitation period is 4400 years. The slow migration with respect 
to the orbital motion affects the eccentricity evolution only through the oscillation period. 
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Fig. 7.— Escape of a companion located near the keplerian boundary of a conservative 
constant-direction acceleration with Iq = 30°. The distances are given in AU. The boundary’s 
semi-major axis is ak p i r = 100 AU. Note how the companion hovers above the star before 
escaping. The corresponding orbital elements are shown in Figure (8). 





36 





Fig. 8.— Escape of a companion located near the keplcrian boundary of a conservative 
constant-direction acceleration. The boundary’s semi-major axis is ak p i r = 100 AU. The 
companion is initially on a circular orbit of semi-major axis 68.5 AU with Iq = 30°. The semi¬ 
major axis experiences large amplitude librations about ctk p ir- The eccentricity amplitude is 
larger than sin Jq and the inclination evolution differs from that inside the keplerian region. 
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Fig. 9.— Migration and eccentricity excitation near the keplerian boundary of a finite dura¬ 
tion acceleration for two inclination values J 0 = 0 (upper panel) and J 0 = 20° (lower panel). 
In each panel, the final semi-major axis (if (AU) (solid) and the final eccentricity ef (dashed) 
are shown as a function of the initial semi-major axis a\ (AU) as given by equations (29) 
and the numerical integration of the full equations of motion (2) (symbols). The parameters 
are: V = 0.35 km s^ 1 , ak p i r = 300 AU, and r = 500 years. The inclined circular orbits were 
started at the descending node (6 = 180°). For J 0 = 0°, a^ = 7341AU and for J 0 = 20°, 
a out = 3 477 AU, a esc = 3 797AU, = 14 542AU. 
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Fig. 10.— Final semi-major axes some of the initially inclined orbits (Jo = 20°) of Figure 
(9) as a function of orbital longitude 0. The dashed curve corresponds an initial semi-major 
axis at a esc . 




Fig. 11.— The evolution of the eccentricity of Jupiter (solid) and Saturn (dashed). In the 
top panel mutual interactions are included but no external acceleration is present. The 
planets’ initial eccentricities are ej = 0.04 and es = 0.07. In the bottom panel, an external 
acceleration with ctkpir = 10 3 AU acts on the two planets initially on circular coplanar orbits 
while mutual interactions are turned off. 
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Fig. 12.— The excitation of Jupiter’s and Saturn’s eccentricities and inclinations by a 
constant-direction acceleration with ak p i r = 10 3 AU including the planets’ mutual gravi¬ 
tational interactions.The planets were started with circular co-planar orbits. From top to 
bottom, the panels show the evolution of the eccentricities, Jupiter’s inclination (solid) and 
the mutual inclination (dashed), Jupiter’s argument of perihelion (solid) and the mutual ar¬ 
gument of perihelion (dashed), Jupiter’s longitude of ascending node (solid) and the mutual 
longitude of ascending node (dashed). 
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Fig. 13.— The excitation of Jupiter’s and Saturn’s eccentricities and inclinations by a con¬ 
stant direction acceleration with ak p i r = 10 4 AU including the planets’ mutual gravitational 
interactions. The planets were started with circular co-planar orbits. From top to bottom, 
the panels show the evolution of the eccentricities, Jupiter’s inclination (solid) and the mu¬ 
tual inclination (dashed), Jupiter’s argument of perihelion (solid) and the mutual argument 
of perihelion (dashed), Jupiter’s longitude of ascending node (solid) and the mutual longitude 
of ascending node (dashed). 
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Fig. 14.— Eccentricity excitation, apsidal alignment and binary migration in the v And 
system. The acceleration is of the form A 2 (dotted) with Aq = 2.8 x 10 -11 kms -2 , r = 
2000 years, t 0 = 12000 years corresponding to ak p i r = 300 AU and a residual velocity V = 
5.6 kms -1 . From the top down, The first panel shows the eccentricity excitation of planets 
d (solid) and c (dashed) and the eccentricity evolution of v And B (short-dashed) as well as 
the acceleration pulse normalized to its maximum value (dotted). The second panel shows 
the radial migration of v And B. The third panel shows the inclination of planets d (solid, 
left scale) and c (dashed, left scale) and the inclination evolution of v And B (short-dashed, 
right scale). The last panel shows the relative apsidal libration of planets d and c. 


















